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Abstract — Non-Euclidean geometry combined with transfor- 
mation optics has recently led to the proposal of an invisibility 
cloak that avoids optical singularities and therefore can work, 
in principle, in a broad band of the spectrum [U. Leonhardt 
and T. Tyc, Science 323, 110 (2009)]. Such a cloak is perfect in 
the limit of geometrical optics, but not in wave optics. Here we 
analyze, both analytically and numerically, full wave propagation 
in non-Euclidean cloaking. We show that the cloaking device 
performs remarkably well even in a regime beyond geometrical 
optics where the device is comparable in size with the wavelength. 
In particular, the cloak is nearly perfect for a spectrum of 
frequencies that are related to spherical harmonics. We also show 
that for increasing wavenumber the device works increasingly 
better, approaching perfect behavior in the limit of geometrical 
optics. 

Index Terms — Electromagnetic cloaking. Transformation op- 
tics. Wave propagation 

I. Introduction 

ELECTROMAGNETIC cloaking is a modern and rapidly 
developing area of optics. The main objective of cloaking 
is to encapsulate an object into a cloak, a carefully designed 
material with special optical properties, such that the object 
becomes invisible. This is achieved by bending light rays 
(or wave fronts) such that they flow smoothly around the 
object and after passing it, restore their original direction, 
creating the illusion that the space is empty. Several types 
of cloaks have been proposed 0,El,|[3l,||l,|[3,|[6]|,||7l, 
some of which have been realized experimentally (H, ||9l, 
Cni, imi. All are based on transformation optics III, O, O, 
El, lEl, the idea that optical media effectively 
transforms the geometry for light, and all, except Q, require 
metamaterials |16|, |17|, |18|, |19|, |20| for their practical 
realization. Originally, two different approaches have been 
used to propose an invisibility cloak. 

The proposal by J.B. Pendry fTl used a geometric trans- 
formation of Euclidean space to expand a single point into 
a region of finite volume. Based on ideas of transformation 
optics, a metamaterial was designed that would lead light 
around this region, making it invisible. A drawback of this 
proposal was the necessity of infinite light speed at the 
boundary of the invisible region, which prevents such a device 
from working in a broad band of the spectrum. The proposal 
was realized experimentally for microwaves fSl and provided 
an inspiration for other similar ideas. 

The proposal by U. Leonhardt O also applied transfor- 
mation optics, but used methods of complex analysis for 
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designing an invisibility cloak that could be made with an 
isotropic medium. In particular, it took advantage of multi- 
valued analytic functions and the rich structure of Riemann 
sheets and branch cuts in the complex plane. An important 
ingredient was a refractive index distribution with focusing 
properties such that a ray that left the main Riemann sheet 
could return back to it and continue in the original direction. 
However, the speed of light in this proposal diverged as well, 
so it could not work broadband either. 

The paradigm of cloaking was changed by the proposal of 
U. Leonhardt and T. Tyc |3|. They showed that the necessity 
of infinite light speed can be avoided by further developing the 
ideas of | 2 1 and proposed a cloaking device based on a combi- 
nation of transformation optics and non-Euclidean geometry. 
All optical parameters of the material have non- singular values 
in the whole space and therefore nothing prevents them, in 
principle, from working broadband. This proposal may turn 
invisibility into a practical technology. In contrast to Euclidean 
cloaking 1 1 1, however, non-Euclidean cloaking works perfectly 
only in the limit of geometrical optics. For wave optics, 
imperfections emerge due to diffraction and phase delays. 
These imperfections do not matter, though, in most practical 
situations, because they become negligible if the size of the 
device is much larger than the wavelength. At the same time, 
it would be very useful to know what the character of the 
imperfections is and how the cloaking device behaves if the 
validity condition for geometrical optics is not satisfied. This 
is the purpose of this paper. We present an analysis of wave 
propagation in the 2D version of the cloaking device ||3 
and provide computer simulations of wave propagation in the 
cloak. 

The paper is organized as follows. In Sec JIl] we explain 
the idea of transformation optics and in Sec. [llT| we describe 
the non-Euclidean cloaking device. In Sec.llVlwe calculate the 
optical properties of the material of the cloak and in Sec. [V] we 
analyze theoretically wave propagation in the cloak. In Sec|Vl| 
we support our theoretical analysis by numerical simulations 
of the wave propagation in the cloak and we conclude in 
Sec.lVlIl 

II. Transformation optics 

In this section we briefly and intuitively explain the idea 
of transformation optics, referring to the papers 1131, ifHIl 
for the full technical details and techniques of calculating the 
material properties of the optical medium. For simplicity, we 
will consider the limit of geometrical optics and ignore the 
polarization of light. 

Transformation optics relies on the fact that an optical 
medium effectively changes the geometry of space perceived 
by light. What determines light propagation is not the geo- 
metric path length, but rather the optical path length, which 
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should be, according to Fermat's principle, stationary (usually 
minimal, but in some cases it is maximal) for the actual path 
of the ray. The optical path element is n-times larger than 
the geometric path element, where n is the refractive index. 
Hence, the medium creates an effective "optical metric tensor" 
for light that equals -times the "geometric metric tensor", 
and light rays then propagate along geodesies with respect to 
this optical metric tensor. In case of an anisotropic refractive 
index (such as in crystals or metamaterials), the transformation 
of the metric is just slightly more complicated. 

Conversely, suppose that light propagates in an empty space 
(so-called virtual space) R = (X, F, Z) with a uniform 
refractive index where ray trajectories are straight lines. If we 
map this space by some function / to a space r = {x^y^z) 
(so-called physical space) such that r = f{R), the geometric 
metric is not preserved in general, because dx^ + dy'^ + dz'^ ^ 
dX^ + dF^ + dZ^ . However, filling physical space with a 
suitable optical medium, we can compensate for this and 
make the optical metrics in virtual and physical spaces equal. 
Then in physical space light will follow the images /(/) 
of straight lines / of virtual space. By choosing a suitable 
transformation /, one can achieve various interesting effects 
such as invisibility |1| or perfect focusing to a single point 
of light rays coming from all directions [21 J. In general the 
mapping / is not conformal (i.e., it does not isotropically 
stretch or squeeze an infinitesimal neighborhood of a given 
point) and therefore the resulting refractive index in physical 
space is not isotropic. 

Remarkably, transformation optics works much more gener- 
ally than just in the simplified situation we have just described. 
In particular, it is valid not just in the limit of geometrical 
optics, but it works equally well for full electromagnetic waves 
(HI. If the electric permittivity and magnetic permeability are 
appropriately chosen 1141 . Maxwell's equations in physical 
space are equivalent to Maxwell's equations in empty virtual 
space. It turns out |14| that the tensors of relative electric 
permittivity and relative magnetic permeability must be pro- 
portional to each other to avoid birefringence and therefore to 
ensure that both light polarizations behave the same way in 
the medium. Ideally, the two tensors are equal, which gives 
impedance matching with the vacuum | 22 | and hence causes 
no reflection at the interfaces of the device. 

III. Non-Euclidean cloaking device 

Next we describe the non-Euclidean cloaking device O 
using pictures and try to give some intuition for the geometric 
transformation used. We will leave the mathematical details 
for the next section; a full description of the transformation 
including all calculations can be found in |23|, the Supporting 
Online Material of |3|. 

The most distinctive feature of the cloaking device |3 | is 
that virtual space is not just the simple Euclidean space as 
in other proposals but contains a non-Euclidean part. We will 
focus on the two-dimensional version of the cloak here; the 
generalization to three dimensions is relatively straightforward 
and has been described in 13|, 1231 . 

Virtual space consists of two parts: (i) a Euclidean plane P 
and (ii) a 2-sphere S (surface of a 3D ball). The two parts are 



connected along a branch cut, a ID line that on the sphere 
corresponds to one quarter of its equator while in the plane 
it corresponds to a straight line segment. The best way to 
visualize this is to imagine that the plane is partly wrapped 
around the sphere so that P and S touch each other along / 
(see Fig. [T]). 




Fig. 1. Virtual space consists of a sphere and a plane that are connected 
along a line segment - branch cut. 

A light ray propagating in the plane P that hits the branch 
cut (at the point X, say) transfers to the sphere S. There 
it propagates along a geodesic, which is a great circle, and 
returns to the point X where it transfers back to the plane P 
and continues in the original direction. In this way, the sphere 
is invisible and the only feature that can possibly disclose its 
presence is a time delay of the rays that have entered the 
sphere. 




Fig. 2. When a mirror (gray thick line) is placed along a great circle (equator 
in the picture), any light ray (yellow line) still returns to its original position 
at the branch cut (black line) where it entered the sphere but never visits the 
northern hemisphere; it therefore becomes invisible. 

To create an invisible region where an object could be 
hidden, one can take advantage of the high symmetry of 
the sphere. Placing a mirror (see Fig. [2]) along a suitably 
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chosen great circle on the sphere, the Hght ray, after two 
reflections, stiU returns to the point X but completely avoids 
the hemisphere hidden behind the mirror; anything placed 
there becomes invisible. 




Fig. 3. Physical space is a plane; the sphere S of virtual space is mapped 
to the inside of the red curve while the plane P of virtual space is mapped 
to the outside of the curve. The hidden hemisphere is mapped to the inside 
of the gray curve which becomes the invisible region of the device. 

The physical space of the device is a Euclidean plane filled 
with an optical medium. The map from virtual to physical 
space is based on bipolar coordinates and maps both the 
sphere S and the plane P of virtual space to regions of the 
physical plane (see Fig.[3|. The map contains no singularities, 
i.e., no points where space would be infinitely stretched or 
squeezed. Therefore the speed of light in the medium is finite 
and nonzero everywhere, which is a great advantage compared 
to Euclidean cloaks. The invisible region is the image of the 
hidden hemisphere of virtual space. 

IV. Material properties of the cloak 

In this paper we focus on a two-dimensional version of 
the non-Euclidean invisibility cloak that has been described 
in the previous section. To adjust the two-dimensional situ- 
ation to real propagation of electromagnetic waves that are 
governed by Maxwell equations and are therefore intrinsically 
three-dimensional, we do the following construction. We per- 
form the above described geometric transformation from two- 
dimensional virtual space to the plane xy of physical space 
and then we add an extra dimension z, making physical space 
three-dimensional. In calculating the line elements in this 3D 
physical space, we then have to add the term dz^ to the 
line elements corresponding to the two-dimensional situation. 
We also assume that the wavevector of the incoming wave is 
perpendicular to the z-axis and hence the phase and amplitude 
of the wave do not depend on z. This makes the problem 
effectively two-dimensional, just as required. 

To describe the Euclidean part of the cloak, we use bipolar 
cylindrical coordinates (r, a^z) in physical space that are 
related to the Cartesian (x^y^z) as 

a sinhr a sincr 

y=-—rz — 77—^ z = z. (1) 



Here the parameter a determines the size of the cloak. Anal- 
ogous relations hold in virtual space where x^y and a are 
replaced by x\y' and a' , respectively. The mapping between 
the two spaces is given in the bipolar cylindrical coordinates 
by the following relation between a and cj''. 



for |cr| < 



G = 



3|cr| + TT I sgn cr for 



(2) 



(3) 



and the coordinates r, z coincide in the two spaces. The line 
element ds^ = dx'^ + d.y'^ + d.z^ in virtual space is equal to 

^2 



ds^ = (dr^ + Ag'') + Az\ (4) 

(cosnr — cos g'Y 

The permittivity and permeability tensors for this region can 
be then calculated using the same general method 1.14 J as in 
ES with the result 



inj = diag 



da 
d^ 



dG^ 

da' 



(cosh r — cos a) da' 



(5) 



(coshr — cos a') da 

The mapping between the non-Euclidean part of virtual 
space and the corresponding part of physical space is more 
involved. We will not reproduce here all the corresponding 
formulas of f23l, but rather mention just the most important 
steps. One starts with a sphere of radius 

4 

r = — a (6) 

TT 

(see Eq. (S43) of | 23 1) parameterized by spherical coordinates 
- the latitudal angle and the longitudinal angle that is denoted 
by a^ This is not yet the sphere of of virtual space, but is 
related to it by a suitable Mobius transformation (Eq. (S40) 
of 1231 ) represented on the sphere via stereographic projection 
(Eqs. (S37)-(S39) of |23|). Then the coordinates (a^6>) are 
mapped to the bipolar coordinates (a, r) of physical space as 
follows. The mapping between a' and a is given by Eq. ^ 
where g' is set to one of the intervals [— 27r, — tt] and [tt, 27r] 
by adding or subtracting 27r. The corresponding intervals of a 
are then [— tt, —Stt/A] and [37r/4,7r]. The mapping between 
and r can be derived by combining Eqs. (S41) and (S46) of 
|23|. This yields after some calculations 



(9 = 2 arctan(t - 1 + y^f^ + 1) 



with 



t = tan 



sinhi 



4 V cosh r + 1 



7r/4 



(7) 



(8) 



The line element of the non-Euclidean part of virtual space is 



d^^ 



l + cot2((9/2) 



2 cos a' cot((9/2) 
X (di9^ + si 



2cot2(i9/2)_ 
^(9da'^)+d^^ (9) 



cosh r — cos a ' 



cosh r — cos a 



This line element is slightly different from the line element 
(S47) in |23|. In fact, the latter turned out to be wrong, because 
it is the line element on the sphere {X^Y'Z^) as described in 
section "Non-Euclidean branch" of 1.23 J , while the correct line 
element should be measured on the sphere {XYZ). The factor 
in the rectangular parentheses in Eq. ^ expresses the scaling 
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of the line element when going from one sphere to the other, 
as can be calculated from the transformation (S40) of 1231 
between the two spheres. 

After some algebra, Eq. ^ can be brought into the form 



1 + (t - 1 + v^Tiy 



(coshr + 1)^ 



2(t - 1 + Vt^Ti) 



da'^ + dz^ (10) 



with 



1 + (t - 1 + 

- 1 + ^/W^f + 1 



(t - 1 + yjt^ + 1)2 - 2(t - 1 + yjt^ + 1) cosa' + 2 ' 

(11) 

Following again the general method |[T4ll . we finally obtain 
the permittivity and permeability tensors for the region 37r/4 < 
\g\ < tt in bipolar cylindrical coordinates: 



with 



TT + 1 + tVt^Tl 



da 



4 t-l + VPTI coshr + ld(j'' 
1 



16 (t - 1 + v^^Ti)(t2 + 1 + tVWTi) 

^ [l + (t-l + v/t^)2]2 

(cosh r — cos cr) 2 da' 
cosh r + 1 dcr 



(12) 



(13) 



(14) 



(15) 



Figure |4] shows the values of jicj , /i^- and Sz in the cloaking 
region, respectively. We picked these material parameters, 
because in our simulations we used transverse electric (TE) 
waves whose electric field is parallel with the z-axis, and 
in this case the material parameters that influence the wave 
propagation are only /i^r, Mr and Sz- We found numerically that 
^cj e [0.190, 1], iir e [1,5.274] and Sz e [1, 179.15]. Clearly, 
all the values avoid zero and infinity and hence non-Euclidean 
cloaking requires no optical singularity of the material. 

To create the invisible region, a mirror is placed in virtual 
space along a great circle on the sphere. Since an ideal mirror 
can be described by a material with infinite permittivity, it 
retains this property upon a geometric transformation and 
remains an ideal mirror. Therefore in physical space we simply 
place the mirror along the curve that is the image of the great 
circle. 

There is also another option for creating an invisible region 
apart from using the mirror, namely expanding a line of 
physical space that is not crossed by any ray. We will not 
consider this option in this paper. 

V. Wave propagation in the invisibility cloak 

First we will focus on the situation where there is no 
invisibility region yet; there is no mirror and the whole sphere 
of virtual space is accessible to light. Then we proceed to the 
situation with the mirror. 




Fig. 4. The values of /io- (a), /ir (b) and Ez (c) in the region of the 
cloaking device; the black line marks the border between the Euclidean and 
non-Euclidean region. All the values are non- singular. 



To see what happens when a light wave propagates in the 
cloak, we need to understand the behavior of the waves in 
the non-Euclidean region of the device. This region is an 
image of a sphere in virtual space with uniform refractive 
index. Wave propagation on such a sphere is different from 
that in an Euclidean space, and the solutions of the wave 
(or, equivalently, Helmholtz) equation are described by the 
spherical harmonics Yim{0^ (p). The situation here is still more 
complicated, though. The sphere in virtual space is not closed 
but is "cut" along the branch cut, and the two sides of the 
branch cut are connected to different parts of the plane of 
virtual space. Therefore the amplitude and phase of the wave 
on the two sides of the branch cut need not be the same as 
would be the case of a normal, connected sphere. The result is 
that there can be also other waves propagating on this sphere 
beside the spherical harmonics. 

It can happen, though, that the incoming wave in the plane 
is tuned in such a way that that the wave on the sphere will 
naturally have the same amplitude and phase at any two points 
adjacent across the branch cut, so it will behave as a wave 
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on a normal, uncut sphere. Evidently, this can only happen 
for waves that are resonant on the sphere, i.e., for the modes 
described by spherical harmonics. 

To find the corresponding frequencies, we recall that the 
spherical harmonics Yim{0^ (p) are solutions of the Helmholtz 
equation Au-\-k'^u = on the sphere. Here / is a non-negative 
integer and m is an integer satisfying — / < m < I. The 
corresponding quantized wavevectors satisfy the relation 



v7(ITi) 



(16) 



where r is the radius of the sphere, and depend on / but not 
on m. This degeneracy with respect to m is related to the high 
symmetry of the sphere. The corresponding frequencies and 
wavelengths are then 



27rr 



(17) 



We therefore expect that the cloak will behave perfectly 
for light waves with wavevectors and frequencies given by 
Eqs. ([16]) and ([17]). 

Another way of seeing this is the following. A spherical 
harmonic describes a wave running on the sphere. If we 
make a loop on the sphere and return to the same place, 
the value of Yim returns to its original value. In this sense 
we can say that a wave described by a spherical harmonic 
interferes constructively with the version of itself that has 
traveled around the sphere. And this is exactly what happens 
to a wave that has entered the sphere via the branch cut - 
it travels on the sphere, is delayed, and then leaves via the 
other side of the branch cut for the plane again. If such a 
wave should interfere constructively with the wave that has 
not entered the sphere at all, the time delay must be equal to 
a multiple of the period of the wave. Since this happens to 



waves with wavevectors and frequencies given by Eqs. (16) 
and ( [Tt] ), these will be the waves for which the cloak will 
work perfectly. 

For waves whose wavelengths do not satisfy the condi- 
tion ( [Tt] ), the interference of the waves that entered the 
sphere with those that did not will not be constructive and 
therefore we can expect that the device will not work perfectly, 
disturbing the waves. 

Note that the constructive interference condition ( [TT] ) is 
different from what one might expect by the following (wrong) 
argument: The waves that enter the sphere will interfere with 
those that did not and the result of this interference will depend 
on the time delay corresponding to propagation around the 
sphere. At first sight, it might seem that this time delay is 
simply equal to the length of the equator of the sphere 27rr 
divided by speed of light c. This would suggest that the device 
will perform best if this delay is an integer multiple of the 
period of the light wave 27r/cj, which would then yield the 
wavelength of light equal to A = 27rr/m with m integer. 
However, this is different from the correct condition ( [Tt] ). It 
is the curvature of the sphere that modifies wave propagation 
compared to a flat, Euclidean space and makes such a simple 
argument invalid. 
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Fig. 5. Wave propagation in the cloaking device for / = 1 (a), / = 2 (b), 
I = 3 with various angles of incidence (c, d, e), and for / = 4 (f). Due to 
constructive interference the waves are not disturbed. 



Next we consider the full version of the non-Euclidean 
cloak including the mirror. In the limit of geometrical optics, 
the mirror should make no difference, because a ray entering 
the sphere of virtual space at some point returns to the same 
point just as if no mirror were present. Moreover, even from 
the semi-geometrical optics point of view there should be no 
change, as the ray picks up a total phase of 27r: one tt for each 
of the two reflections from the mirror. 

However, looking at the situation with the mirror fully wave- 
optically, there should be some change caused by the mirror. 
Indeed, the mirror imposes a boundary condition on the wave, 
in particular, the electric field must be zero at the mirror 
surface. A more detailed analysis shows that instead of 2/ + 1 
linearly independent solutions Yim{0^(p) of the Helmholtz 
equation for a given there are just / independent solutions 
that vanish on a given great circle. For example, if the great 
circle is chosen as the equator (corresponding to 6 = 7r/2), 
then these solutions will be the Yim themselves, but just those 
with m having different parity than 1. Therefore the waves 
propagating in the mirror cloak do not have such a freedom 
as they would if no mirror were present, because of the 
boundary condition or, equivalently, because of the restricted 
set of modes. We can then expect that the mirror cloak will 
not work perfectly even for those wavelengths satisfying the 
condition ( [TT] ). 

VI. Numerical simulations 

Now we proceed to numerical simulations of wave prop- 
agation in the cloaking device to see whether they confirm 
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Fig. 6. Wave propagation at the cloaking device with the condition flT) not 
satisfied for integer I. The values of / are 2.8 (a), 3.2 (b), 3.5 (c, d, e) with 
various angled of incidence and 1.5 (f). The device does not work perfectly 
due to destructive (or partially destructive) interference of the waves that did 
and that did not enter the non-Euclidean part of the device. We see that for 
half-integer values of I the phase shift of the waves that entered the non- 
Euclidean part of the device is approximately tt, which is obvious especially 
in the pictures (d) and (e). 



the above considerations. Figure |5] shows the propagation 
of a wave incident on the cloak with wavelength satisfying 
Eq. ([T7| with several different values of / and several angles 
of impact. We see that the wave is not disturbed by the cloak 
and propagates as if no cloak were present at all. 

Furthermore, Fig. [6] shows the propagation of waves that 
do not satisfy Eq. ( p7| ) with integer /. To compare the two 
cases, we still use, however, the number / as in Eq. ([17]), 



but we allow also for its non-integer values. We see from 
Fig. |6] that for non-integer / the waves are disturbed, just as 
we predicted. The perturbation is strongest for half-integer 
values of /, because they are most distant from the resonant 
integer values. Figure 6 clearly shows the time lag of the 
waves propagating through the cloaking structure, resulting in 
phase shifts that are strongest for half-integer I. Nevertheless, 
the emerging phase fronts are still parallel to the incoming 
fronts, which indicates that images are faithfully transmitted. 
However, due to the time lag, two wavefront dislocations 
are formed, as Fig. 6 also shows. They would cause image 
distortions at the edges of the non-Euclidean region of the 
device. 

Next we consider the complete cloaking device with mirror. 
Figure [7] shows the behavior of waves with integer I. Although 
there is some disturbance of the waves, the perturbation is 



Fig. 7. Same as Fig. [5] that is, / = 1 (a), / = 2 (b), / = 3 (c, d, e), 
/ = 4 (f) but with the mirror included. The mirror is the ellipse-like object 
that encapsulates the invisible region. Clearly, electric field vanishes in this 
region. Due to the restricted set of modes the performance of the device is 
no more perfect, but the imperfection is weak. 



quite modest. This is quite remarkable if we take into account 
that the non-Euclidean cloak has originally been proposed for 
the regime of geometrical optics, and there was absolutely no 
ambition that it works for wavelengths comparable with the 
size of the device. Here we see that the cloak can perform very 
well also outside of geometrical optics, and large wavelengths 
do not seem to be a problem. Furthermore, Fig. 8 shows that 
the mirror does not affect the propagation of waves with non- 
integer values of /. The wave distortions are very similar to 
the case without the mirror. 

Last we look at the behavior for large values of I. Un- 
fortunately, the simulations became soon very difficult and 
even on a powerful computer converged just up to / = 20. 
We therefore could not numerically test the behavior of the 
device for larger values. Figure [9] shows the case / = 20 and 
/ = 19.5 without and with the mirror. We see that while for 
/ = 20 the device works nearly perfectly (we think that the 
imperfections are caused by numerical errors), for / = 19.5 
the behavior is almost perfect, apart from interference effects 
on the edges. Extrapolating this for larger we can conclude 
that the cloaking device should indeed work increasingly well. 

vn. Conclusion 

In conclusion, we have analyzed wave propagation in the 
two-dimensional version of the non-Euclidean cloaking de- 
vice |3 1. We have seen that although the device works perfectly 
in the geometrical optics limit, there are imperfections due 
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Fig. 8. Same as Fig.[6] that is, / = 2.8 (a), I = 3.2 (b), / = 3.5 (c, d, e) and 
/ = 1.5 (f), but with the mirror included. The restricted set of modes causes a 
slightly different performance of the device compared to the situation without 
the mirror. 



to interference effects if full wave description is employed. 
In the situation where there is no hidden region, there exist 
wavelengths of light for which the device is absolutely perfect. 
This results from constructive interference between the waves 
that entered the non-Euclidean part of the cloak and those that 
did not. This constructive interference occurs for wavelengths 
that resonate on the non-Euclidean part of the cloak and 
correspond to frequencies of spherical harmonics. For other 
wavelengths the interference is not fully constructive and the 
waves leaving the device are disturbed — the device is not 
perfect for light waves. The perturbations are caused by the 
time lag in the device and result in wavefront dislocations. 

Although adding the mirror and creating an invisible region 
causes the cloak to cease working perfectly for the resonant 
wavelengths, the imperfection is very weak. The performance 
of the device is much better than originally expected O. When 
the wavelength gets smaller, the performance of the cloaking 
device gradually improves, becoming perfect in the limit of 
geometrical optics. 

The fact that the non-Euclidean invisibility cloak does not 
work perfectly for light waves is by far outweighed by the 
absence of optical singularities and the potential of work- 
ing broadband. Although non-Euclidean cloaking has been 
designed for wavelengths much smaller than the size of the 
device, which is also the case of most practical situations, 
it is useful to understand the behavior of the device outside 
this limit, as we have tried to do in this paper. This might 
help designing other types of non-Euclidean cloaks that could 
possibly work even for wavelengths comparable with the size 
of the device. 
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Fig. 9. Wave propagation in the cloak without the mirror (first line) and 
with it (second line) for I = 20 (a), (c) and / = 19.5 (b), (d). The device 
now performs quite well even for the half-integer 1. 



References 

[I] J. B. Pendry, D. Schurig, and D. R. Smith, "Controlling Electromagnetic 
Fields", Science, vol. 312, pp. 1780-1782, May 2006. 

[2] U. Leonhardt, "Optical Conformal Mapping", Science, vol. 312, pp. 
1777-1780, May 2006. 

[3] U. Leonhardt and T. Tyc, "Broadband Invisibility by Non-Euclidean 
Cloaking", Science, vol. 323, pp. 110-112, Jan. 2009. 

[4] H. Y. Chen and C. T. Chan, "Transformation media that rotate electro- 
magnetic fields", Appl. Phys. Lett, vol. 90, p. 241105, Jun. 2007. 

[5] Y. Lai, H. Y. Chen, Z.-Q. Zhang and C. T. Chan, "Complementary Media 
Invisibility Cloak that Cloaks Objects at a Distance Outside the Cloaking 
Shell", Phys. Rev. Lett., vol. 102, p. 093901, Mar. 2009. 

[6] J. Li and J. B. Pendry, "Hiding under the Carpet: A New Strategy for 
Cloaking", Phys. Rev. Lett., vol. 101, p. 203901, Nov. 2008. 

[7] I. I. Smolyaninov, V. N. Smolyaninova, A. V. Kildishev, and V. M. 
Shalaev, "Anisotropic Metamaterials Emulated by Tapered Waveguides: 
AppHcation to Optical Cloaking", Phys. Rev. Lett., vol. 102, p. 213901, 
2009. 

[8] D. Schurig, J. J. Mock, B. J. Justice, S. A. Cummer, J. B. Pendry, A. 
F., and D. R. Smith, "Metamaterial Electromagnetic Cloak at Microwave 
Frequencies", Science, vol. 314, pp. 977-980, 2006. 

[9] R. Liu, C. Ji, J. J. Mock, J. Y. Chin, T. J. Cui, D. R. Smith, "Broadband 
Ground-Plane Cloak", Science, vol. 323, pp. 366-369, 2009. 

[10] J. Valentine, J. Li, T. Zentgraf, G. Bartal, X. Zhang. (2009, Apr.). Di- 
electric Optical Cloak. [Online]. Available: http://arxiv.org/abs/0904.3602 

[II] L. H. GabrieUi, J. Cardenas, C. B. Poitras, M. Lipson, (2009, 
Apr.). Title: Cloaking at O ptical Frequencies. [Online]. Available: 
|http://arxiv.org/abs/0904.3508| 



8 



[12] L. S. Dolin, "On the possibility of comparing threedimensional electro- 
magnetic systems with nonuniform anisotropic fillings", Isv. Vusov, vol. 
4, pp. 964-967, 1961. 

[13] A. Greenleaf, M. Lassas and G. Uhlmann, "On nonuniqueness for 
Calderon's inverse problem". Math. Res. Lett., vol. 10, pp. 685-693, Sep. 
2003. 

[14] U. Leonhardt and T. G. Philbin, "General relativity in electrical engi- 
neering". New J. Phys., vol. 8, p. 247, 2006. 

[15] U. Leonhardt and T. G. Philbin, "Transformation Optics and the Geom- 
etry of Light", Prog. Opt., vol. 53, 2009. 

[16] V. M. Shalaev, "Transforming fight". Science, vol. 322, pp. 384-386, 
2008. 

[17] G. W. Milton, The Theory of Composites. Cambridge University Press, 
Cambridge, 2002. 

[18] D. R. Smith, J. B. Pendry, M. C. K. Wiltshire, "Metamaterials and 

negative refractive index". Science, vol. 305, pp. 788-792, Aug. 2004. 
[19] C. M. Soukoulis, S. Linden, M. Wegener, "Negative refractive index at 

optical wavelengths". Science, vol. 315, pp. 47-49, Jan. 2007. 
[20] A. K. Sarychev, V. M. Shalaev, Electrodynamics of Metamaterials. 

World Scientific, Singapore, 2007. 
[21] U. Leonhardt and T. Tyc, "Superantenna made of transformation media". 

New J. Phys., vol. 10, p. 115026 (2008). 
[22] J. D. Jackson, Classical Electrodynamics. Wiley, 1999. 
[23] U. Leonhardt and T. Tyc. (2009, January). Broadband Invisibility by 

Non-Euclidean Cloaking. Science. Supporting online material of |3|. 

[Online]. Available: http ://w w w. sciencemag . org/cgi/data/ 116633 2/DC 1 / 1 ^ 




Che Ting Chan received his BS degree in Physics 
from the University of Hong Kong in 1980 and 
his PhD degree from UC Berkeley in 1985. He is 
currently a Chair Professor of Physics at HKUST, 
and the director of William Mong Institute of Nano 
Science and Technology (WMINST). Dr. Chan's 
research interest has been in the area of photonic 
crystals, metamaterials and nano-materials. Dr. Chan 
is a Fellow of the American Physical Society. 



Tomas Tyc was born in Brno, Czechoslovakia (now 

i Czech Republic) on February 7th, 1973. He re- 
ceived his Master degree and Ph.D. in theoretical 
physics from Masaryk University in Brno in 1996 
and 1999, respectively. In his research Tomas Tyc 
initially investigated correlations in free electron 
beams, later he switched to quantum optics and 
quantum information theory. In 2007 he focused on 
optics, in particular on theory of invisibile cloaks 
and transformation optics. He obtained habilitation 
(Associate professorship) from Masaryk University 
in 2006. He was a Research Fellow at University of Vienna (2000), Macquarie 
University Sydney (2001, 2002), University of Calgary (2004), and University 
of St Andrews (2007, 2008). Beside his research interests, he likes to 
popularize science and makes shows on physics of everyday life for general 
public. 




l^^^^^M Ulf Leonhardt was born in Schlema, in former 

^m^^^^^M East Germany, on October 9th 1965. He studied 
-^W^^^ at Friedrich-Schiller University Jena, Germany, at 

I 1 Moscow State University, Russia, and at Hum- 

f^S5k^J boldt University Berlin, Germany. He received the 

Diploma in Physics from Friedrich-Schiller Univer- 
sity in 1990 and the PhD in Theoretical Physics 
from Humboldt University in 1993. Ulf Leonhardt 
was a research associate at the Max Planck Research 
Group Nonclassical Radiation in Berlin 1994-1995, 
a visiting scholar at the Oregon Center for Optics in 
Eugene, Oregon, 1995-1996, a Habilitation Fellow of the German Research 
Council at the University of Ulm, Germany, 1996-1998, and a Feodor- 
Lynen and Gran-Gustafsson Fellow at the Royal Institute of Technology 
in Stockholm, Sweden, 1998-2000. Since April 2000 he is the Chair in 
Theoretical Physics at the University of St Andrews, Scotland. In 2008 he was 
a Visiting Professor at the National University of Singapore. Ulf Leonhardt 
is the first from former East Germany to win the Otto Hahn Award of the 
Max Planck Society. For his PhD thesis he received the Tiburtius Prize of the 
Senate of Berlin. In 2006 Scientific American listed him among the top 50 
policy business and research leaders for his work on invisibility devices. In 
2008 he received a Royal Society Wolfson Research Merit Award. He is a 
Fellow of the Institute of Physics and of the Royal Society of Edinburgh. 



Huanyang Chen was born in Fujian, China, on 
August 5, 1983. He received his B.Sc. and Ph.D. 
in physics from Shanghai Jiao Tong University, 
Shanghai, China, in 2005 and in 2008. He was a 
Research Assistant and a Post-Doctoral Fellow in 
the group of Professor Che Ting Chan at Department 
of Physics, the Hong Kong University of Science 
and Technology, Hong Kong, China, from 2006 to 
2009. He will join Suzhou University, Jiangsu, China 
after September 2009. His research interests include 
photonic band gaps, negative index materials, meta- 
material design and transformation optics. He has published several original 
papers on transformation optics and cloaking. 




